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This  work  investigates  thermally  induced  stresses  and  deformations  in  a  multilayered  thin  film  thermo¬ 
electric  (TE)  structure  for  assessing  thermomechanical  reliability  of  the  TE  devices.  The  multilayered 
structure  consists  of  a  n-type  and  p-type  thermocouple  separated  by  an  insulating  layer,  and  an 
additional  supporting  membrane  layer.  The  one-dimensional  thermoelectricity  is  used  to  compute  the 
temperature  distribution  with  heat  transfer  by  radiation  and  convection  at  the  film  surfaces  being 
considered.  The  thermoelasticity  theory  of  laminated  composites  is  employed  to  determine  the  stresses 
and  deformations  in  the  thin  film  structure.  It  is  found  that  the  thermal  stress  calculated  using  the  pres¬ 
ent  laminate  model  has  a  significantly  higher  magnitude  than  that  predicted  by  the  strength  of  materials 
model.  Effects  of  heat  transfer  coefficient  at  the  film  surfaces  are  also  examined. 

©  2014  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

Thermoelectric  (TE)  materials  have  been  considered  a 
promising  material  technology  for  converting  waste  heat  to 
electric  power.  The  energy  conversion  efficiency  of  TE  materials 
is  described  by  the  so-called  figure  of  merit  ZT  =  S2T/pk,  where  S 
is  the  Seebeck  coefficient,  p  the  electrical  resistivity,  k  the  thermal 
conductivity  and  T  the  temperature.  Generally  speaking,  conven¬ 
tional  TE  devices  have  macroscopic  dimensions  and  low  figure  of 
merit,  which  renders  them  inapplicable  in  microelectronics  appli¬ 
cations.  In  recent  years,  thin  film  thermoelectric  materials  for 
microelectronics  applications  have  been  developed  with  enhanced 
energy  conversion  efficiency  represented  by  the  increased  figure  of 
merit,  see,  for  example,  Rowe  et  al.  [1],  Volklein  et  al.  [2], 
Venkatasubramanian  et  al.  [3],  Chen  et  al.  [4],  Bottner  et  al.  [5,6], 
and  Takashiri  et  al.  [7], 

Thermoelectrics  are  subjected  to  thermal  gradients  under  oper¬ 
ation  conditions.  The  thermal  gradient  is  produced  by  both  the 
temperature  difference  applied  at  the  hot  and  cold  ends  of  the 
thermoelectric  element  and  the  Joule  heating  due  to  the  thermo¬ 
electric  effects.  The  thermal  gradient  induces  thermal  stresses 
which  may  cause  failure  of  the  thermoelectric  structure.  Assess¬ 
ment  of  thermomechanical  reliability  of  thin  film  thermoelectric 
devices  requires  a  thorough  understanding  of  thermal  stresses  in 
the  structures.  Thermal  stresses  in  thermoelectric  thin  films  have 
been  investigated  by  Huang  et  al.  [8,9],  They  computed  the  thermal 
stresses  in  a  four-layered  thermoelectric  element  using  both  a 
one-dimensional  (ID)  strength  of  materials  model  and  a  two¬ 
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dimensional  (2D)  finite  element  numerical  method,  and  concluded 
that  the  thermal  stresses  of  the  ID  and  2D  models  differ  in  the  end 
regions  of  the  structure.  The  ID  strength  of  materials  model  may 
be  used  to  estimate  the  thermal  stress  in  the  preliminary  design. 
Mahan  [10],  Volklein  et  al.  [2],  and  Zabrocki  et  al.  [11]  also  showed 
that  the  ID  model  is  appropriate  in  analyzing  the  temperature  field 
in  thermoelectrics.  Jin  [12]  calculated  the  temperature  that  causes 
buckling  failure  of  a  thin  film  thermoelectric  structure  using  the 
membrane  force  derived  from  the  strength  of  materials  model. 

Multilayered  thin  film  thermoelectric  structures  may  have 
properties  and  configurations  that  are  asymmetric  about  its  geo¬ 
metrical  mid-plane.  More  accurate  determination  of  the  stresses 
and  deformations  in  the  structure  requires  employment  of  lami¬ 
nated  composites  theories  which  have  not  been  addressed  in  the 
literature  to  the  author’s  best  knowledge.  This  work  presents  an 
analytical  model  to  compute  the  stresses  and  deformations  in  a 
multilayered  thermoelectric  structure  using  a  theory  of  laminated 
composites.  Because  the  length-to-thickness  ratios  of  thin  film 
thermoelectrics  may  reach  an  order  of  1000  [2,13],  the  classical 
laminate  theory  may  be  appropriate  for  evaluating  the  thermal 
stresses  and  deformations  in  the  structure.  The  remainder  of  this 
paper  is  organized  as  follows.  Section  2  presents  the  basic  equa¬ 
tions  of  one-dimensional  thermoelectricity  and  obtains  the  tem¬ 
perature  field  in  a  four-layer  thin  film  structure  consisting  of  a  n- 
type  and  p-type  thermocouple  separated  by  an  insulating  layer, 
and  an  additional  supporting  membrane  layer.  Section  3  first  de¬ 
scribes  the  equations  of  the  classical  laminate  theory  for  the  mul¬ 
tilayered  thermoelectric  structure  and  subsequently  derives  the 
equations  for  computing  the  thermal  stresses  and  deformations. 
In  Section  4,  numerical  results  are  presented  to  illustrate  the  differ¬ 
ences  between  the  predictions  of  the  laminate  model  of  this  work 
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and  the  strength  of  materials  model.  The  effects  of  heat  transfer  at 
the  film  surfaces  on  the  thermal  stresses  are  also  examined.  Finally, 
conclusions  are  provided  in  Section  5. 


2.  Temperature  field 


Consider  a  multilayered  thin-film  thermoelectric  (TE)  structure 
as  shown  in  Fig.  1.  The  structure  is  a  four-layer  TE  element  consist¬ 
ing  of  a  n-type  and  p-type  thermoelectric  couple  (layer  #2  and  #4), 
an  insulating  layer  (layer  #3)  and  a  supporting  membrane  (layer 
#1).  The  structural  configuration  is  the  same  as  that  studied  by 
Huang  et  al.  [8,9]  which  is  an  idealized  model  of  the  thermoelectric 
structure  fabricated  by  Volklein  et  al.  [2],  It  is  assumed  that  the  two 
thermoelectric  films  have  the  same  properties  and  thickness. 
Moreover,  the  insulating  and  supporting  layers  have  the  same 
thickness.  The  geometric  mid-plane  of  the  structure  is  thus  the 
interface  between  the  n-type  and  insulating  layers.  We  note  that 
layers  of  different  thicknesses  can  be  considered  similarly  as  long 
as  the  geometric  mid-plane  of  the  structure  is  located  at  an  inter¬ 
face  between  the  layers.  Furthermore,  if  the  geometric  mid-plane  is 
not  at  an  interface,  the  problem  can  still  be  treated  by  assuming 
the  layer  containing  the  geometric  mid-plane  as  two  layers  of 
the  same  properties  perfectly  bonded  along  the  geometric 
mid-plane. 

Volklein  et  al.  [2]  and  Huang  et  al.  [8,9]  analyzed  the  tempera¬ 
ture  distribution  in  the  four-layer  TE  element  assuming 
one-dimensional  (ID)  steady  state  heat  flow  in  the  longitudinal 
direction  (x-direction).  This  section  first  summarizes  the  basic 
equations  of  ID  thermoelectricity  and  subsequently  revisits  the 
temperature  field  in  the  TE  element  using  boundary  conditions 
appropriate  for  TE  generators.  The  temperatures  at  the  cold  and 
hot  ends  of  the  TE  element  are  denoted  by  T0  and  Th,  respectively. 

Consider  ID  steady  state  electrical  and  heat  conduction  in  the 
longitudinal  direction  (x-direction)  of  the  thermoelectric  structure. 
The  basic  equations  for  the  thermoelectric  films  are  [10,14] 

q=-kpa^(+JST  (1) 


dq 

dx 


(2) 


4/ 

dx 


0 


(3) 


where  T  is  the  temperature,  q  the  heat  flux,  J  the  electric  current 
density,  S  the  Seebeck  coefficient,  kpn  the  thermal  conductivity, 
and  p  the  electrical  resistivity.  We  assume  that  the  n-type  and 
p-type  thermoelectric  films  have  the  same  material  properties. 
For  the  support  and  insulation  layers,  the  term  involving  the  electric 
current  should  be  dropped  in  Eqs.  (1)  and  (2).  Clearly,  the  electric 
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current  density  J  is  a  constant  in  ID  steady  state  electrical  conduc¬ 
tion  according  to  Eq.  (3). 

In  deriving  the  governing  equation  of  the  temperature  field,  we 
adopt  the  approach  of  Volklein  at  el.  [2]  by  using  an  effective  heat 
transfer  coefficient,  t],  to  describe  the  combined  heat  transfer  by 
convection  and  radiation  at  the  film  surfaces,  i.e.,  the  heat  flux  in 
the  outward  normal  direction  of  the  film  surfaces,  qn,  is  propor¬ 
tional  to  the  temperature  difference  between  the  film  temperature 
and  the  ambient  temperature.  Mathematically,  this  condition  is 
expressed  as  follows 

qn  =  ri(T-To )  (4) 


Here  we  assume  that  the  ambient  temperature  is  the  same  as  that 
at  the  cold  junction,  T0. 

By  considering  energy  balance  for  a  small  differential  element 
in  the  longitudinal  direction  of  the  multilayered  structure,  we 
can  get  the  governing  equation  of  the  temperature  in  the  structure 
as  follows 


fd2(T-T0) 

dx2 


£(T-T0)  =  -pJ2 


where  k  is  an  effective  thermal  conductivity  defined  by 


(5) 


(6) 


in  which  kis  and  ksu  are  the  thermal  conductivities  of  the  insulating 
and  supporting  layers,  respectively,  and  h,  (1=1,  2,  3)  the  thick¬ 
nesses  of  the  supporting  membrane,  the  n-type  thermoelectric  film 
and  the  insulating  layer,  respectively.  Note  that  hi  =  h3  as  assumed. 

The  solution  of  Eq.  (5)  under  the  boundary  conditions  at  the 
cold  and  hot  ends  (T  =  T0  at  x  =  0,  and  T  =  Tt,  at  x  =  L)  can  be  obtained 
as  follows 

T  -  T0  =  A  exp(Xx)  +  B  exp(-Xx)  +  ^  pj2  (7) 

where  X  is  a  parameter  given  by 


and  A  and  B  are  constants  given  by 

A  = - —— 1 - — —  (  AT  -  [1  -  exp(-2L)l  —  pj2 ] , 

exp(2L)  -  exp(-2L)  \  1  h  i 

B  = - tt-tt— — — { AT  +  [exp(2L)  -  1]  —  pj2\ 

exp(2L)  -  exp(-U)  \  1  '  h  i 


in  which  L  is  the  length  of  the  layered  structure  and  AT  is  the  tem¬ 
perature  difference  between  the  hot  and  cold  junctions,  i.e., 
AT=Th-T0. 


3.  Thermal  stresses 

The  classical  lamination  theory  of  layered  structures  may  be 
appropriate  for  evaluating  the  in-plane  thermal  stresses  and  defor¬ 
mations  in  thin  film  thermoelectric  structures  because  the  length- 
to-thickness  ratios  of  thin  film  thermoelectrics  may  reach  an  order 
of  1000  [2,13],  In  the  laminate  model,  the  displacements  U,  V  and 
W  in  the  x-,  y-  and  z-directions,  respectively,  are  given  by  [15,16] 

U  =  u(x,y)-z^, 

V  =  v(x.y)  >  (10) 

W  =  w(x,y) 

where  u(x,y)  and  t(x,y)  are  the  displacements  of  the  geometrical 
mid-plane  (z  =  0)  of  the  multilayered  structure.  As  a  first  order 
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approximation,  we  assume  that  the  displacements  vary  only  in  the 
longitudinal  x-direction  as  the  temperature  T  is  only  a  function  of  x. 
This  is  the  so  called  cylindrical  bending  which  can  accurately  cap¬ 
ture  the  normal  thermal  stress  in  the  longitudinal  direction  that 
may  cause  compressive  or  buckling  failure,  one  of  major  concerns 
in  the  design  for  thermomechanical  reliability  of  the  thin  film  struc¬ 
ture.  Another  failure  mode  is  delamination  caused  by  the  shear 
stresses  between  the  neighboring  layers  of  the  TE  structure  [9], 
Accurate  determination  of  the  shear  stresses  requires  employment 
of  general  bending  theory  which  will  be  addressed  in  the  future 
study.  Now  the  basic  equations  for  the  mid-plane  displacements 
u,  v  and  w  (as  functions  of  x  only)  reduce  to  [15,16] 


~dx2 


d(T-T0) 

dx 


djv 

dx2 


0 


Du 


d4w 

dx4 


d3u 

d? 


r  d\T-T0) 
”  dx2 


The  decoupled  form  of  the  above  equations  for  u  and 
obtained  as  follows 


(Ha) 

(lib) 

(lie) 
can  be 


Du  J  dx3  V  ”  Dn  *7  dd 


(12a) 


BjA  <fw 

/  dx 4 


\d2(T-T„) 
)  dx2 


where  An,  Bn,  Du,  and  are  given  by 


(12b) 


(13) 


“  2 1  —  Vfc  k  *  iy 

in  which  Ek,  vk  and  ak  are  Young’s  modulus,  Poisson’s  ratio  and  the 
coefficient  of  thermal  expansion  of  the  kth  layer,  respectively,  and 
zk  are  the  z-coordinates  of  the  plate  surfaces  and  the  interfaces  be¬ 
tween  the  layers  defined  by 


z0  =  -(h1  +h2), 

Z\  =  -h2, 

z2  =  0, 

z3  =  h3  =  h2, 

z4  =  h3  +  h4  =  hi  +  h2 


The  boundary  conditions  for  the  displacement  field  are 


at  x 

=  0, 


=  0,  L 
at  x  =  0,  L 


for  the  clamped  ends,  and 


(15) 


(16) 


u  =  v=  0,  at x  =  0,  L 
w  =  0,  ^  =  0,  at  x  =  0,  L 


(17) 


for  the  pinned  ends. 

It  follows  from  the  basic  Eq.  (lib)  and  the  boundary  conditions 
for  the  displacement  vin  Eqs.  (16)  and  (17)  that  v=0.  Solving  the 
boundary  value  problem  described  by  Eqs.  ((12),  (16),  (17))  for  the 
mid-plane  displacements  u  and  w,  we  have 


+  dw  (19) 


where  T  and  T  are 


T=Jo  {T -T0)dx=j[exp(Ax)  - 1]  +  ^[1  -  exp(-foc)]  +  ^pj2x, 


(20) 


and  the  integration  constants  au,  bu,  aw,  bw,  cw  and  dw  are  deter¬ 
mined  as  follows 


buL2  «  -  (n^  +  LT(L )  -  ■ 3  ^  [2T(L)  -  LT(t)] 


(21a) 

awL3  =  +  [2T(L)  -  LT(L)], 

bwL 2  =  -2  (m^  +  N^)  [3 T(L)  -  LT(L)} 

(21b) 

cw  =  dw  =  0 

(21c) 

for  the  clamped  ends,  and 

au~Dnaw’ 

(22a) 

0w=  +  bw  =  0 

(22b) 

Cw  =  -a^L2,  dw  =  0 

(22c) 

for  the  pinned  ends. 

A  major  concern  in  the  reliability  design  of  thin  film  thermo¬ 
electrics  is  the  normal  stress  in  the  longitudinal  direction.  The 
stress  in  the  kth  layer  (k  =  1,  2,  3,  4)  is  given  by 

T  (& +  Z<)  -  (T  -  To)  (23) 

where  the  normal  strain  ej*  and  the  curvature  k ^  of  the  mid-plane 
are 

(*u  -§^e°**=  fa  ) (T  -  To)  +  aux  +  bu 


(24) 
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|d„  - j  <  =  -  (Mit  +  ^  (T  -  To)  -  awx  -  bw  (25) 

Finally,  the  thermal  stresses  in  the  strength  of  materials  model 
is  given  by  [12] 

0**  =  EkE%"m)  -  Ekak(T  -  T0)  (26a) 

where 

4mm)  =  [|>a  I  E£'hJ  [(T  -  T°)  -  ^r]j|  (26b) 

Eqs.  (23)  and  (26)  show  that  the  normal  stress  prediction  of  the 
laminate  model  differs  from  that  of  the  strength  of  materials  model 
by  considering  the  Poisson’s  ratio  effect,  the  curvature  of  the  de¬ 
formed  thin  film  structure,  and  the  extension-bending  coupling  ef¬ 
fect  on  the  normal  strain.  These  differences  are  illustrated  by  the 
numerical  examples  in  the  following  section. 

4.  Numerical  results  and  discussion 

In  the  numerical  calculations,  we  use  the  material  parameters 
listed  in  Table  1  [2,8],  The  temperatures  at  the  cold  and  hot  ends 
are  T0  =  290  K  and  Th  =  330  K,  respectively.  The  temperature  depen¬ 
dence  of  the  properties  is  not  considered.  All  layers  are  assumed  to 
have  the  same  thickness,  i.e.,  hi  =  h2  =  h3  =  h4  =  1  pm,  and  the 
length  is  assumed  to  be  L  =  300  pm  in  the  numerical  examples. 
Finally,  we  use  the  following  electrical  current  density  that  approx¬ 
imately  corresponds  to  the  peak  conversion  efficiency 

Jpk  = - ^  (J\  +Z[T„  -  (AT/2)(kpn/k)]  -  l) 

p  SL[Th-(AT/2)(kpn/k)]V  P  J 

(27) 

where  Z  =  S2/(pkpn)  is  the  dimensional  figure  of  merit. 

Fig.  2  shows  the  thermal  stress  a x*  in  the  longitudinal  direction 
of  the  thin  film  structure  with  clamped  ends.  The  Young’s  modulus 
of  the  insulating  layer  equals  75  GPa  and  the  heat  transfer  at  the 
film  surfaces  is  not  considered,  i.e.,  rj  =  0.  The  thermal  stresses  in 
the  p-type  and  n-type  thermoelectric  layers  are  about  the  same  be¬ 
cause  the  two  films  have  the  same  material  properties  and  the  cur¬ 
vature  effect  is  not  significant  (see  Fig.  7)  as  indicated  by  Eq.  (23).  It 
is  seen  from  the  figure  that  the  thermal  stresses  are  compressive  in 
the  support  membrane  layer  and  the  thermoelectric  layer.  The 
compressive  stress  in  the  support  membrane  increases  from  the 
cold  end  to  its  maximum  of  about  31  MPa  at  the  hot  end.  The  ther¬ 
mal  stress  in  the  n-type  film,  however,  remains  approximately  flat 
over  the  length.  Finally,  the  thermal  stress  in  the  insulating  layer  is 
compressive  near  the  cold  end  and  becomes  tensile  near  the  hot 
end. 

Fig.  3  shows  the  thermal  stress  a **  in  the  multilayered  structure 
using  the  strength  of  materials  model  as  adopted  in  the  existing 
studies  [8,12],  All  the  material  and  geometrical  properties  are  the 


Table  1 

Material  properties. 


Thermoelectric  films 
(n-type  and  p-type) 

Insulating 

Support 

Young’s  modulus  (GPa) 

160 

25,  75 

138 

Poisson’s  ratio 

0.30 

0.22 

0.32 

Coefficient  of  thermal 

4.2 

0.5 

5.9 

expansion  (106/K) 

Thermal  conductivity 

2.0 

17.3 

(W/(m  K)) 

Electrical  resistivity  (O  m) 

0.6  x  10-5 

Seebeck  coefficient  (pV/K) 

200 
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Fig.  2.  Thermal  stresses  in  the  longitudinal  direction  of  the  multilayered 
structure  with  clamped  ends  (Young’s  modulus  of  the  insulating  layer  £3  =  75  GPa, 
the  effective  heat  transfer  coefficient  at  the  film  surfaces  r/  =  0). 


Fig.  3.  Thermal  stresses  cr**  in  the  longitudinal  direction  of  the  multilayered 
structure:  strength  of  materials  model  (Young’s  modulus  of  the  insulating  layer 
£3  =  75  GPa,  the  effective  heat  transfer  coefficient  at  the  film  surfaces  r\  =  0). 


same  as  those  used  in  Fig.  2.  The  overall  characteristics  of  the  stress 
field  are  similar  to  those  shown  in  Fig.  2.  The  magnitude  of  the 
stress  in  each  layer,  however,  is  much  smaller.  For  example,  the 
magnitude  of  the  maximum  compressive  stress  in  the  membrane 
is  now  21  MPa  which  is  about  2/3  of  that  using  the  more  advanced 
laminate  model  shown  in  Fig.  2.  Therefore,  design  of  multilayered 
thin  film  thermoelectric  structures  for  thermomechanical  reliabil¬ 
ity  should  be  based  on  the  laminate  model  instead  of  the  strength 
of  materials  model. 

Fig.  4  shows  the  thermal  stresses  <rxx  in  the  multilayered  struc¬ 
ture  with  clamped  ends.  The  Young’s  modulus  of  the  insulating 
layer  is  now  reduced  to  25  GPa.  All  other  material  and  geometrical 
properties  are  the  same  as  those  used  in  Fig.  2.  Compared  with  the 
thermal  stresses  in  the  structure  containing  an  insulating  layer 
with  a  modulus  of  75  GPa  shown  in  Fig.  2,  the  thermal  stress  is  re¬ 
duced  significantly  in  the  insulating  layer  while  deviates  a  little  bit 
from  those  in  other  layers  in  Fig.  2.  The  implication  is  that  the  insu¬ 
lating  layer  may  be  designed  as  a  non-load-bearing  layer  without 
significantly  affecting  the  mechanical  integrity  of  the  overall 
structure. 

The  effect  of  heat  transfer  at  the  film  surfaces  on  the  thermal 
stresses  in  the  thermoelectric  structure  is  examined  in  Fig.  5  where 
only  the  stress  in  the  supporting  membrane  is  shown  and  the 
effective  heat  transfer  coefficient  at  the  film  surfaces  is  assumed 
to  be  i]  =  20  J/(m2  s  K).  It  appears  that  consideration  of  surface  heat 
transfer  has  an  insignificant  influence  on  the  thermal  stress  partic¬ 
ularly  because  the  t]  value  used  in  the  calculation  is  probably  larger 
than  those  encountered  in  practical  designs. 
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Fig.  5.  Thermal  stress  ct**  in  the  longitudinal  direction  of  the  membrane  support 
layer  in  the  structure  with  clamped  ends  (Young’s  modulus  of  the  insulating  layer 
£3  =  75  GPa,  the  effective  heat  transfer  coefficient  at  the  film  surfaces  t)  =  0,  20). 


To  ascertain  the  effects  of  extension-bending  coupling  on  the 
normal  stress,  we  further  examine  the  normal  strain  and  curvature 
of  the  mid-plane  of  the  multilayered  structure  with  pinned  ends 
because  pinned  ends  result  in  a  slightly  larger  deflection  (out  of 
plane  displacement  w).  Fig.  6  compares  the  mid-plane  normal 
strain  of  the  multilayered  structure  with  the  corresponding  nor¬ 
mal  strain  in  the  strength  of  materials  model.  All  the  material  and 
geometrical  properties  are  the  same  as  those  used  in  Fig.  4.  It  is 
seen  that  the  maximum  magnitude  of  the  normal  strain  ej,  of  the 
laminate  theory  is  about  25%  higher  than  that  of  the  strength  of 
materials  model. 


Fig.  6.  Mid-plane  normal  strain  eJ,  of  the  multilayered  structure  with  pinned  ends 
and  the  corresponding  normal  strain  in  the  strength  of  materials  model  (Young’s 
modulus  of  the  insulating  layer  £3  =  25  GPa,  the  effective  heat  transfer  coefficient  at 
the  film  surfaces  rj  =  0). 


Fig.  7  shows  the  mid-plane  curvature  k JJ,  of  the  multilayered 
structure  with  pinned  ends.  All  the  material  and  geometrical  prop¬ 
erties  are  the  same  as  those  used  in  Fig.  6.  The  curvature  increases 
from  zero  at  the  ends  of  the  film  and  reaches  the  peak  at  the  mid¬ 
span  of  the  structure.  Although  pronounced  curvature  is  observed, 
its  contribution  to  the  normal  stress  is  insignificant  compared  with 
that  of  the  mid-plane  normal  strain  because  the  coordinate  z  is  on 
the  order  of  a  few  microns,  which  shows  that  the  extension-bend¬ 
ing  coupling  influences  the  normal  stress  mostly  through  its  effect 
on  the  normal  strain. 

5.  Concluding  remarks 

Assessment  of  thermomechanical  reliability  of  thin  film  ther¬ 
moelectric  devices  requires  a  thorough  understanding  of  thermal 
stresses  in  the  devices  because  the  thin  film  structures  are  sub¬ 
jected  to  sever  thermal  gradient  under  operation  conditions. 
Knowledge  of  thermal  stresses  is  also  critical  in  predicting  reliabil¬ 
ity  and  preventing  failure  of  general  microelectronic  devices  [17], 
It  has  been  known  that  compared  with  finite  elements  and  other 
numerical  methods,  the  strength  of  materials  model  is  straightfor¬ 
ward  to  implement  and  still  provides  reasonable  results  for  design 
considerations.  The  laminate  model  presented  in  this  work  not 
only  improves  the  predictions  of  the  strength  of  materials  ap¬ 
proach  but  also  remains  as  a  convenient  and  feasible  tool  to  eval¬ 
uate  thermal  stresses  and  hence  the  structural  integrity  of  thin  film 
thermoelectric  devices.  The  higher  compressive  normal  stress 
compared  with  the  strength  of  materials  prediction  will  likely  in¬ 
crease  the  probability  of  compressive  failure  of  the  semi-conductor 
thin-film  structures.  Moreover,  buckling  failure  may  also  become 
possible  for  thin  films  with  high  length  to  thickness  ratios  under 
the  increased  compressive  thermal  stress. 

We  note  that  although  the  effects  of  extension-bending  cou¬ 
pling  and  Poisson’s  ratio  on  the  thermal  stress  are  included  in 
the  present  model,  the  cylindrical  bending  assumption  does  not 
consider  variations  of  the  field  variables  in  the  y-direction.  Hence, 
a  refined  model  that  accounts  for  the  y-dependence  of  the  field 
variables  should  be  employed  in  the  future  study.  Furthermore,  be¬ 
cause  the  classical  lamination  theory  assumes  zero  transverse 
shear  strains,  a  higher  order  lamination  theory  (e.g.,  [18,19]) 
should  be  used  to  accurately  determine  the  transverse  shear 
strains  and  stresses  which  may  cause  delamination  failure  of  mul¬ 
tilayered  thermoelectric  structures. 
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